Cooling Torsional Nanomechanical Vibration by Spin-Orbit Interactions 
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We propose and study a spin-orbit interaction based mechanism to actively cool down the torsional 
vibration of a nanomechanical resonator made by semiconductor materials. We show that the spin- 
orbit interactions of electrons can induce a coherent coupling between the electron spins and the 
torsional modes of nanomechanical vibration. This coherent coupling leads to an active cooling for 
the torsional modes via the dynamical thermalization of the resonator and the spin ensemble. 

PACS numbers: 85.85.+j, 42.50.Pq, 71.70.Ej, 32.80.Pj 



Introduction.-FdiSt developments of nano-fabrication 
technology enable us to manufacture the nanomechan- 
ical resonators (NAMRs) with high frequency and high 
quality [1]. These artificial systems have attracted more 
and more attentions for their application potentials, such 
as ultra-sensitive probe for tiny displacement [2], single 
spin detector [3], and physical implementations of quan- 
tum information [4]. On the other hand, it is believed 
that a high frequency NAMR (e.g. with GHz oscillation) 
will exhibit various quantum effects in this macroscopic 
system at low temperatures [1]. Therefore, the studies on 
NAMR are relevant to fundamental problems like quan- 
tum measurement [5] . 

To give prominence to the quantization effects of 
NAMR, the crucial issue depends on whether we can cool 
it down to the vibrational ground state. Actually, in- 
teracting with the surrounding environment, the NAMR 
reaches a thermal equilibrium through the relaxation. 
The mean thermal occupation number of the NAMR, 
which is determined by the environment temperature and 
the vibration frequency, is always larger than unity for 
most of the present nanomechanical systems even at the 
sub-Kelvin temperature. Accordingly, besides decreasing 
the environment temperature, special cooling mechanism 
is needed to be invented to further reduce the thermal oc- 
cupation number. 

Recent experiments have shown that the micromechan- 
ical systems can be cooled down by radiation pressure [6] 
or by coupling to a single electron transistor [7]. And sev- 
eral theoretical scenarios were proposed based on differ- 
ent cooling mechanism, such as laser excitation of phonon 
sideband of a quantum dot [8], and periodic coupling to 
a Cooper pair box[9], et al. All the cooling mechanisms 
mentioned above are only focused on the flexural mechan- 
ical modes. In this Letter, we propose an active cooling 
mechanism for the quantized torsional nanomechanical 
modes (TNMMs) [10-12] of the NAMR. 

Similar to the displacements and momentums for the 
flexural modes, the TNMMs are described by the oscil- 
lation of the torsion angles and the angular momentums 
of the NAMR. Thus, it is quite natural to to consider 



how to couple the TNMMs to the electron spin degrees 
of freedom, which also behave as angular momentums. 
Recently, a novel interaction between the spin current 
in semiconductor materials and the NAMR induced by 
the spin-orbit coupling is discovered in Ref. [13]. This 
pioneer investigation is carried out in the semi-classical 
regime to provide the possibility for detecting spin cur- 
rent by the NAMR[13, 14]. 

Here, we present a microscopic quantum mechanical 
description for the coherent coupling between the elec- 
tron spins and the TNMM. With this coherent coupling, 
we show that the mean occupation number of the NAMR 
can be significantly reduced. Thus, the NAMR is actively 
cooled down and can be brought to the quantum realm. 

Model for spin-strain coupling. -We consider a doubly 
clamped nanomechanical rod with length width and 
hight h (Fig.l). Similar to the structure considered in 
Ref. [13], the semiconductor layer grows on an insulat- 
ing substrate (Fig. 1(b)). The electrons flowing through 
the nanomechanical rod are confined in the semiconduc- 
tor layer. For the nth normal mode, the torsion angle 
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FIG. 1: (a) Schematic illustration of torsional vibration of the 
nano- mechanical rod. (b) The cross section of the nanome- 
chanical rod and its strain distribution. Electrons flowing 
through the nanomechanical rod are confined in the semicon- 
ductor layer (blue), which grows on an insulating substrate 
(grey), (c) A typical strain distribution of ezx (red) and tzy 
(blue). It is shown that tzx is an even function with respect 
to X, while Czy is odd. 
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of the NAMR is On{z^t) = OnQsm{n7rz/l) cos{(jJrit) with 
frequency ujri and amphtude Ono- 

According to the isotropic elastic theory [15, 16], the 
strain tensor field has non-vanishing components Cxz = 
ezx = e'{z){df/dy) and e^, = e,y = -0' {z){df / dx) 
where 0'{z) = dO/dz^ and the function f{x,y) is deter- 
mined by the cross section geometry of the NAMR[15]. 
We illustrate the strain field explicitly in Fig. 1(b). 

To describe the quantum phenomena of the tor- 
sional oscillation of the NAMR, we apply the stan- 
dard canonical quantization procedure by modeling the 
NAMR as the many-harmonic-oscillator with Hamilto- 
nian i^NAMR = J2n^n{(^li^n + 1/2). Here, the boson 
operator = I^^'^^/^Auo^j2h{0nQ + iOno/^n) is intro- 
duced through the effective mass M = p\fll/2^ where p 
is the mass density, and / = f{x'^^y'^)dA is the rotating 
inertia about z axis for A being the cross section area. 

In the following, we will focus on the fundamental 
(n = 1) TNMM at low temperatures and thus omit 
the subscript index of the boson operators. Then the 
quantized strain field along the NAMR is described by 
6.. = {df/dy)F{z),e,y = -{df/dx)F{z), where F{z) = 
ei cos{qz){a^-\-a), ei = 7r/-^(;i/2McJi V?) ^2 and q = tt/L 

As pointed out in Ref. [13], the spins in a nanome- 
chanical rod can be coupled to the strain field in narrow 
band semiconductor materials by the Hamiltonian[17] 

i^lo = a[(7x{exzkz - ^xyky) + c.p.], (1) 

where a is the coupling strength, and crx,y,z and hk^^y^z 
are the Pauli matrices and the momentum components 
of the electron respectively. Here, 'c.p.' stands for cyclic 
permutation oix^ y and z. Due to the confinement in the 
X and y directions, the spacial wave function of the elec- 
tron can be separated as ip{x,y,z) = (p{x,y) ex.p{ikzz), 
where (p{x,y) is determined by the transverse confine- 
ment in the x and y directions. We further assume that 
this confinement is so strong that we can take an average 
over (p{x,y) to obtain an effective Hamiltonian Hso = 
{(f{x, y)\H^Q\(p{x^ y)). Considering the parity of the non- 
vanishing strain components (see Fig. 1(c)), the Hamil- 
tonian -f^lo is then reduced from three-dimension to one- 
dimension effectively, i.e. i^so = Oi£kz{a^ +a) cos{qz)ax^ 
where S = {df/dy)ei is a dimensionless constant, and 
(df/dy) denotes the average over the wave function 
ip{x,y). 

Next, we further simplify the single electron Hamil- 
tonian by considering the classical limit of the spatial 
motion of the electron. In the case of large longitudinal 
momentum of the injected electron compared with the 
wave number of the TNMM, i.e. kz » q, the longi- 
tudinal motion of the injected electron is hardly affected 
by the back action of the vibration. Thus, we neglect the 
recoil effect and used the approximate Galilean transfor- 
mation z ^ zq -\- vt^ where v = hkz/m is the velocity of 
the coherent injected electron. With this approximation. 



the Hamiltonian can be written as: 

H= ^/kj;^cr;^ + /kJi a^a + ^^ cos(cc;^t + 0) (a^ + a) cr^, (2) 

where hujz is the Zeeman energy of the electron spin. 
We have ignored the constant term of kinetic energy 
h^k1/2m. Here, we define the velocity dependent cou- 
pling constant = a£mv/h, the frequency oOy = qv, 
and the phase = qzo. 

There are three different characteristic frequencies in 
this system: the electron Zeeman frequency cj^, the 
NAMR vibration frequency cji, and the the spatial mo- 
tion induced frequency ujy. In the interaction picture, 
the Hamiltonian will involve coupling terms with dif- 
ferent frequencies. In the spirit of rotating wave ap- 
proximation, only the terms with the lowest frequency 
are retained. Therefore, we can adjust the injecting ve- 
locity and/or the Zeeman energy of the electrons, so 
that different types of interaction can be obtained[18]. 
For example, if cjy ~ uUz — ^^i, only two terms are re- 
tained, and we get a JC-type Hamiltonian (see Fig. 2(b)), 
Hjc = 9vCi^(^- exp[i(u;^t + 0)]/2 + H.c. 

Having obtained the actively controllable single elec- 
tron Hamiltonian, we will study the interactions between 
the coherently injected ensembles of electron spins and 
the TNMM of the NAMR. Let TV spin polarized elec- 
trons with velocity v be injected into the NAMR, and 
meanwhile, N spin polarized electrons in the opposite 
direction with the velocity —v. In principle, this model 
can be realized by joining the NAMR with two spin po- 
larized electronic source, for example, two ferromagnetic 
leads(Fig.2(a)). Then the Hamiltonian for this system is 

i^sys = hhjWa + ^hjz{Jz + Qz) 

+ gy{t)a\j--Q-) + H.c., (3) 
where we define the collective spin operators [19] 

N N 

1=1 1=1 

N N 

Q. = ^rW, Q± = ^e^^*-r«. (4) 

i=l i=l 

Here, the Pauli matrices cr^*j_ and r^*^ denote the elec- 
tron spins with the velocity v and —v respectively, and 
9v{t) = exp(ia;^t)/2 is the time dependent coupling 
strength. 

Through the Holstein-Primakoff transformation Jz = 
-N + 26t6, J_ = y^b and Qz = -N ^ 2c^c, Q- = 
\/]Vc, the collective excitation of the ensembles of elec- 
tron spins can be characterized by the operators b{c) and 
b\c^), which satisfy the boson commutative relation in 
the large N-low excitation limit. Then the interaction be- 
tween the collective spin excitation and NAMR is mod- 
eled by an interacting two-mode boson system. In the 
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FIG. 2: (a) Polarized electron spins are injected with veloci- 
ties V and —V into the NAMR from two ferromagnetic (FM) 
leads, (b) Controllable coupling: the collective excitation en- 
ergy uoz of electron spins is linearly dependent on the external 
magnetic field 5, and the frequency uov is determined by the 
electron velocity. The red point denotes the resonant condi- 
tion ujy = ujz — ooi are satisfied, (c) Schematic illustration 
of the cooling mechanism. The system contains two interact- 
ing boson modes: the fundamental TNMM (red, left) with 
frequency cji and the collective spin excitation mode (blue, 
right) with frequency a;^. 

interaction picture, 

Hi = hn{aUe-'^' + H.c), (5) 

where d = {h — c)/ \pi is the boson operator associated 
with the collective spin excitations, hVt — ^^N/2gy is the 
Rabi frequency, and /S. = uOz — ooi — uOy is the detuning. 

Cooling mechanism.-Nex.t, we study the active cooling 
mechanism based on the dynamical thermalization of the 
coupled system of the TNMM and ensembles of electron 
spins (Fig.2(c)). To describe a practical situation, we 
assume the TNMM and the ensembles of electron spins to 
interact with two separated thermal reservoirs. The time 
evolution of the system operators a and d are governed 
by the Heisenberg-Langevin equantions[20] 

d = -iQde-'^^ ~ f ^ + 

d = -iQae'^^ ~ 2 ^ ^ ^^^^^ ^^"^ 

where Fa{t) and Fd{t) are noise operators due to the 
reservoirs. We assume the two reservoirs are statistically 
independent, i.e. {F^{t)Fd{t^)) = 0. The correlations 
beteen Fa{t) and F^{t) 

{Fl{t)Fa{t')) = nna5{t-t'), 
{Fa{t)Fl{t'))=^{na^mt-t') (8) 

describe the quantum fluctuations with n being the 
damping constant and fia = [exp{huJi/kBT) — being 
the thermal occupation number of TNMM. The noise 
operators Fd{t) and F^{t) follow similar relations with 



damping constant 7 and thermal occupation number 
nd=[eM^z/kBT)-l]-\ 

It follows from Eqs.(6) and (7) that the motions of 
bilinear quantities A = {a^a-\-aa^), V = {d^d-\-dd^), and 
C = {ad'^ + d^a) are 

A = inC- inC^ - K^A^K. {2na + 1) , (9) 
V = -inC + il^C* - 7P + 7 {2nd + 1) , (10) 

C = iQA - iVtV - i(A^ + 7)C, (11) 

where C = Cexp(zAt). The above set of equations de- 
termines the time evolution of the occupation numbers 
(Fig. 3). At steady state, the mean occupation number of 
the TNMM is 

{a'^a)ss = na-j{fia- fid) C{k, A), (12) 

where £(/>:, 7, A) = 4n'^{i<i + j)/[K.j[{i<i + j)'^ {AQ'^ / + 
1) + A^]] is a Lorentz line shape function. 

The first term in Eq. (12) represents the initial thermal 
occupation number of the TNMM in equilibrium with the 
thermal reservoir without injection electron spins. And 
the second term, which is proportional to the coupling 
strength 1]^, is due to the interaction with the spin po- 
larized electrons. It is clear that, if > ^d, the mean 
occupation number (a^a)ss will be reduced to a value 
lower than n^, or in other words, the TNMM is cooled 
down. Furthermore, a large cooling efficiency demands 
the condition >> Ud- It is equivalence to require that 

(i) the electron Zeeman energy co^ is much larger than the 
TNMM frequency cji, i.e. cOz >> ^i- This is the first 
necessary condition for a high cooling efficiency. Under 
this condition, the mean occupation number at resonance 
(A = 0) can be rewritten as the weighted average of the 
thermal occupation numbers and Ud'- 

{a^a)ss = (1 - /1/2) fia + fif2nd (13) 

with the two ratios fi = 4(^^/(/^7 + 4(^^) and /2 = 
7/(7 + /t:). Thus, we find that the mean occupation num- 
ber {a^a)ss has a neglectably small lower bound fid, and 
a significant reduction of {a^a)^^ with respect to the ini- 
tial thermal equilibrium value needs the following two 
additional requirements: (ii) the coupling between the 
two modes a and d are strong enough compared with the 
dampings, i.e. AQ^ » and (iii) the decay rate of the 
collective electron spin excitations is much larger than 
the one of the TNMM, i.e. 7 >> 

To explain the physical mechanism beneath the above 
protocol, we suppose that the two boson modes are ini- 
tially in equilibrium with their reservoirs of the same tem- 
perature T. Under the condition (i), the initial occupa- 
tion numbers fia » fid (see Fig. 2(c)). The condition 

(ii) ensures the efficient population transfer between the 
two boson modes. For t > 0, the occupation number of 
the TNMM (collective excitation of electron spins) will 
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FIG. 3: Time evolution of the occupation numbers (a'^a)(t) 
and {d^d){t). The occupation numbers of the two boson 
modes reach steady state values when t ^ +00. Inset: the 
steady state occupation number (a^a)ss as a function of the 
detuning A. Parameters used in the calculation are described 
in text. 

decrease (increase) due to the coupling (Fig. 3). Si- 
multaneously, the reservoirs tend to maintain the ther- 
mal occupation numbers fla and fid. In other words, the 
occupation number of the TNMM (collective excitation 
electron spins) will gain from (decay to) the reservoir 
with the rate hi (7). The condition (iii) defines two dif- 
ferent time scales of the thermalization processes. The 
much faster decay rate 7 than guarantees that the net 
effect of the dynamical thermalization is reduction of the 
occupation number. 

Here, we emphasize the importance of the time depen- 
dence of coupling gy{t){see Eq. 3) induced by the spatial 
motion of spins. Generally speaking, the frequency mis- 
match of two interacting boson modes will prohibit any 
efficient couplings. In other words, the requirement (i) 
would block off the occupation transfer between the two 
boson modes if their coupling were time independent. 
Fortunately, in our model, the spatial motion provides 
a third frequency cOy to compensate the large frequency 
difference between coi and 00^- Thus, the generic reso- 
nant condition co^ = ^1 is modified to cOy =00^—^1- In 
this way, the spatial motion ensures an efficient coupling 
under the frequency mismatch condition (i). 

Finally, we discuss the experimental feasibility of our 
cooling mechanism. For a typical NAMR with / = 2/im, 
d 100 nm, and h 80 nm, and the TNMM with 
frequency coi = 300 MHz, the dimensionless constant 
E ^ 10~^ according to its definitions. With the pa- 
rameters of GaAs material, mass density p = 5.3g/ cm^, 
effective electron mass m = 0.067me, the spin-orbit in- 
teraction constant a/h = Ax 10^ m / s, and the velocity of 
injected electrons vf = 1.6 x 10^ m / s, which corresponds 
to Fermi energy Ef = 5 meV, we estimate the coupling 
strength Q = 7 MHz for = 400. The decay rate 7 for 
the collective excitation of electron spins is determined 
by the electron spin coherence time Ti ~ 50 ns, which im- 
plies 7 - 1/Ti = 20 MHz. For the quality factor Q = 10^ 



of the NAMR , and thus Hi = coi/Q = 3 kHz, we find the 
occupation number could reduced from = 130 without 
coupling down to {a^a)ss = 0.06 << 1 at a sub-Kelvin 
temperature T = 300 mK. 

Condusion.-We investigate a spin-TNMM coupling 
model for cooling the TNMM through dynamical ther- 
malization. We found that the strain induced spin-orbit 
interactions for conduction band electrons in semiconduc- 
tor materials result in a microscopic coherent coupling 
between the electron spins and the quantized TNMM. 
With this discovery, we present an experimentally ac- 
cessible cooling method for the TNMM by injecting en- 
sembles of electron spins. The TNMM can be sufficiently 
cooled down to be capable of exhibiting various quantum 
phenomena. 

This work is supported by the NSFC with grant 
Nos. 90203018, 10474104, 10574077, and 60433050, 
and NFRPC with Nos. 2006CB921206, 2006CB0L0601, 
2006AA06Z104 and 2005CB724508. 



* Electronic address: zjl-dinp@tsinghua.edu.cn 
^ Electronic address: suncp@itp.ac.cn, http://www.itp. 
ac . cn/~suncp 

[1] K. C. Schwab and M. L. Roukes, Physics Today 58, 36, 

(2005) . 

[2] R. G. Knobel and A. N. Cleland, Nature 424, 291, 
(2003). 

[3] D. Rugar, et a/., Nature 430, 329, (2004). 

[4] J. Eisert,et a/., Phys. Rev. Lett. 93, 190402, (2004). 

[5] V. B. Braginsky, Quantum Measurement (Cambridge 

University Press, Cambridge, 1992). 
[6] S. Gigan et a/., Nature 444, 67, (2006). O. Arcizet 

et al, Nature 444, 71, (2006). D. Kleckner and D. 

Bouwmeester, Nature 444, 75, (2006). 
[7] M. D. LaHaye, et al, Science 304, 74, (2004). A. Naik et 

a/., Nature 443, 193, (2006). 
[8] I. Wilson-Rae, P. Zoller and A. Imamoglu, Phys. Rev. 

Lett. 92, 075507, (2004). 
[9] P. Zhang, Y. D. Wang and C. P. Sun, Phys. Rev. Lett. 

95, 097204, (2005). 
[10] L Tittonen et al, Phys. Rev. A, 59, 1038, (1999). 
[11] P. Mohanty et al, Phys. Rev. B, 66, 085416, (2002). 
[12] T. Cohen-Karni et aL, Nature Nanotechnology 1, 36, 

(2006) . 

[13] A. G. Mal'shukov, et al, Phys. Rev. Lett. 95, 107203, 
(2005). 

[14] P. Mohanty, al, Phys. Rev. B 70, 195301, (2004). 

[15] L. D. Landau and E. M. Lifshitz , Theory of elasticity 
(Pergamon Press, New York, 1986). 

[16] A. N. Cleland, Foundations of Nanomechanics (Spinger- 
Verlag, Berlin, 2003). 

[17] G. E. Pikus and A. N. Titkov, in Optical orientation , 
edited by F. Meier and B. P. Zakharchenya (Elsevier Sci- 
ence Pub. Co , Amsterdam , 1984), p. 73. 

[18] Y.-X. Liu et a/., Phys. Rev. Lett. 96, 067003 (2006). 

[19] Z. Song et a/., Phys. Rev. B 71, 205314 (2005). 

[20] M. O. Scully and M. S. Zubairy, Quantum Optics (Cam- 
bridge Univercity Press, New York, 1997). 



